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Abstract 

We consider the effect of radiative corrections on the maximum velocity of propagation of neutral 
scalar fields in a uniform electromagnetic field. The propagator of neutral scalar fields interacting 
with charged fields depends on the electromagnetic field through charged particle loops. The kinetic 
terms of the scalar fields are corrected and the maximum velocity of the scalar particle becomes 
greater or less than unity. We show that the maximum velocity becomes greater than unity in a 
simple example, a neutral scalar field coupled with two charged Dirac fields by Yukawa interaction. 
The maximum velocity depends on the frame of reference and causality is not violated. We discuss 
the possibility of this superluminal propagation in the Standard Model. 
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I. INTRODUCTION 



Classical field equations (Klein-Gordon equation, Maxwell equation, etc.) prohibit su- 
perluminal propagation. However, as shown in Maxwell equation in curved spacetime 
is modified by radiative corrections and the photons travel at speeds greater than unity in 
some background gravitational field without any inconsistency. In this paper we consider 
the effect of radiative corrections on the maximum velocity of propagation of neutral fields 
in a uniform external electromagnetic field in Minkowski spacetime. We show that radiative 
corrections in a uniform electromagnetic field modify the kinetic terms of the field equation 
and the maximum velocities become greater or less than unity depending on the structure 
of the quantum field theory even in Minkowski spacetime. Moreover we show that the maxi- 
mum velocity becomes greater than unity in a simple example, a neutral scalar field coupled 
with two charged Dirac fields by Yukawa interaction. 

In order to make the effect of radiative corrections as clear as possible, we consider 
propagation of a neutral real scalar field interacting with charged fields. The scalar field 
depends on the electromagnetic field only through charged particle loops. 

First we consider generally the effect of radiative corrections and derive the electromag- 
netically modified field equation for the neutral scalar field. We need only consider the two 
points correlation function because we consider only the propagation of individual scalar 
particles. From the translational invariance and the gauge invariance, the connected part of 
the propagator of the scalar field is given by, 

G{x -y,F)^ {Q, F\ T0(x)0(y) \Q, F)^ = / , . . , (1) 



(27r)4p2_^2 + n(p2) + n(p,F)^ 

where \fl,F) is the vacuum state in a uniform electromagnetic field, F^i, = const, and 
{n,Fm{x)(P{y)\n,F)^ = {Q,Fm{x)(P{y)\Q,F) - {Q,F\^{x)\n,F){n,F\(P{y)\Q,F) . 
n(p2) is the self energy for F^^, = and Il{p, F = 0) = 0. From the gauge invariance 
G{x — y) depends on F^^ rather than A^. In general, {fl, F\ \Q, F) 7^ when F^,^ 7^ 0, 
but it has no effect on the propagation and we can omit it. From the Lorentz covariance we 
obtain, 

G{A{x-y),F^^)=G{x-y,F), (2) 

where A is a Lorentz transformation, i.e. {A{x — y))^ = A^^{x — yY and {F^^Y^ = A^^A'^^F^"' . 
From ([1]) and ([2]), Il(p, F) is a scalar function constructed by and F^^^. The general form 



of n(p, F) is given by 

n(p, F) = {aj + a2h)p^ + bj + b2h + ciF^p^p, + 0{F\p^), (3) 
where / = = = E^-B^, h = -iF^^F^, = 2E-B, F^, = ^e^.p^F"-, and 

n 

, ^ s ~ 

{F'^Y^ = F^^^F^^^ ■ ■ ■ F^^'^. n(p, F) in ([3]) is the most general form because of the identities 
{FFY" = ^7]^'^ and [F'^Y^ = (^F'^Y^ ~~ fv^"- For simphcity, we assume that the mass (yu) of 
the scalar field is much smaller than the mass scale (M) of the charged fields (/i <^ M) and 
we restrict within \p^\ <^ M. The higher order terms of p^ in n(p, F) appear as {p^/MY 
and we neglect them. And we assume \eF^^\/M'^ ^ 1 and retain the lowest order terms of 
F^j_^. (The components F^^j depend on the coordinate system (i.e. inertial frame), however, 
when |/i|/M^, \f\/M^ ^ 1 there are coordinate systems in which \eF^y\/M'^ ^1.) 

We consider the shift of the pole at p"^ = fi"^. From the renormalization condition 11 (p^ = 
/i^) = and ^n(p^)|p2=^2 = 0, so we can neglect n(p^) around ~ /x^. Then, when 

^ /i^ we obtain 

p'-fi' + n(/) + U{p, F)^p'-fi' + [a J + a2h)p' + bj + b^h + c^iF^Y^P^P. 

(4) 

~ (1 + aj + a2h){p^ + ciF^Y'Pf^P^) - (/^' - bif - b^h). 

The dispersion relation is given by setting (jlj) equal to zero. From (jlj) we see that the 
modified equation of motion for the scalar field is given by 

[(1 + aj + a2h){d' + ciF^Y'd^d,) + fi"]<P{x) = 0, (5) 

where /i'^ = — bif — 62^- From the kinetic term in ([5]) we see that the spacetime metric 
for the scalar field changes as 

^ (1 + aj + a2hYri^'' + ci{FY1 = (1 + ^i/ + a2h)S^'\ (6) 

Notice that {1 + aif + a2h)S^^ is not universal as the spacetime metric because ai,a2 and Ci 
depend on the coupling constant of the interaction between the scalar field and the charged 
fields. (1 + ai/ + a2h)S^'^ is considered as the spacetime metric only for the scalar field. 

Next we consider the maximum velocity of the scalar particle, i.e. the group velocity 
[vf = dE/dpi) for large momentum {\pi\ ^ /i). We can neglect the mass term when we 
consider the maximum velocity. The kinetic term in ([5]) determines the maximum velocity. 
We transform F^^ into standard forms by Lorentz transformation. There are two cases as 

3 



follows: (i) E and B are parallel to each other, (ii) E and B are mutually perpendicular and 
|E| = |B|. First we consider the case (i). In this case we assume E along the 1-axis and 
obtain = — (F^)^^ = E^, (F^)^^ = (F^)^^ = B^ and the other components are zero. 

We substitute this F^ into the kinetic terms in and obtain 

S^^d.d, = (1 + c^E'){d', - dl) - (1 - ciB'){di + dl). (7) 

From ([7]) we obtain the dispersion relation po = {pi + j^^^^ipl +^3))^^^- Here we neglect 
the mass. Then, the maximum velocities along the parallel direction to E (fy) and along 
the perpendicular direction to E {v±) are given by, 

^--^--(^)"-'-|(^^ + «^'- («' 

In the case (ii) we assume E along the 2-axis and B along the 3-axis and E2 = B3 = |E|. 
We obtain (F^)°° = (F^)°^ = (F^)^° = (F^)^^ = E^ and the other components are zero and 

S'^'-d.d, = (do + + c,E^)do - (1 - ciE2)9i) - d^, - dl (9) 

From (Q we obtain the maximum velocities along the positive and negative 1-axis (f 1+ and 
t>i_) and along the 2, 3-axis {v2 and fs), 

l-ciE2 „ / 1 ci , , , 

^--^'^-^TT^-^-^^^^'^^-^^nTT^J -i-^E^. (10) 

From ([8]) and ffTOj) . in both cases the maximum velocities become greater (less) than unity 
when Ci is negative (positive). We will calculate Ci in an simple example and show that Ci 
can be negative. 

Notice that causality is not violated even if the maximum velocity is greater than unity. 
The crucial point is that the maximum velocity depends on the components of F^ and 
on the frame of reference. There is the assumption of relativistic invariance of the laws 
govennng the mot.on i„ the usual _t about the vtolation of causaltty by fastet thau 
light transmission of information [l!]. In our case this assumption does not hold because 
the equation of motion depends on the components of F^. We show the region to which a 
observer A a.t x = can send a signal in the case (i) in Fig. [l](a). In another inertial frame 
x' this region extends to t' < (Fig. [H^b)). However, from Fig. [11(b) we see that a return 
signal back to A never arrives in the past light cone of A and causality is not violated. In 
the same way causality is not violated also in the case (ii). 




FIG. 1. (a) The region to which an observer ^ at x = can send a signal in the case (i). The 
dotted Unes are hght cones, (b) The same region in another inertial frame x' . A return signal back 
to A never arrives in the past light cone of A and causality is not violated. 

II. A SIMPLE EXAMPLE 

As a simple example, we consider a neutral real scalar field coupled with two charged Dirac 
fields by Yukawa interaction. (Our notation follows Peskin and Schroeder j^.) We calculate 
explicitly c\ in ([3]) in this example and show that ci can be negative. The Lagrangian is 
given by 

where and ?/'2 have a same charge e and Z^^ = 5^ + i&A^J_. The Lagrangian is gauge in- 
variant and invariant under the transformation — — 0, ^/^i — > —1^1,1^2 ~^ So the terms 
(p, (f)^, (pipiipi and 4''ip2'^2 are prohibited by the symmetry and this Lagrangian is renormal- 
izable. We assume g^, e^, A ^ 1 and use the perturbation theory. We use the dimensional 
regularization. The d-dimensional propagators of charged fermions satisfy the equation, 

{ip-m)SA{x,x',m) = 6'^{x-x'), (12) 

where m is either mi or m2. For a uniform electromagnetic field, the exact solution of ( 1T2|1 
is obtained by proper time method pj] (see for review [3,]). The result is 

Sa{x,x' ,m) = ^{x,x'){~i)[i^x + rn]g{x — x',m), (13) 

where 

$(x,x') =exp[-ze / d^^'A^iO] (14) 

Jx' 
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FIG. 2. The lowest order self energy diagram which depends on Ffj_^. The double line stands for 
the propagator {p + m)g{p,m) in p!7|) . 



and 

g{x,m) = / ds- -ttt — 7777 exp [ — tr ln[(eFs)~^ sinh(eFs)] xeF coth{eFs)x 

Jo (47rz)'*/^ s"^/^ 2 4 ^^^^ 

- ^eF^^af^yS - i{m^ - ie)s\, 

where a^y = |[7^,7i/]. The path of ^ integral in the phase factor $( straight 
hne from x' to x. The role of is to make Sa{x,x') gauge covariant, i.e. when 

A^{x) —7- Af^{x) — ^dfj,a{x), ^{x,x') — j- e*("*^'^)~°*^^'))$(x, x'). In f[T^ we use following matrix 
notations; trF" = {F'^Y^ and xF^x = x^{F^'Y''Xy. 



For loop calculations, it is convenient to use the momentum representation of S"^ {4], 

SA{x,x',m) = ^x,x'){-z) I ^e-'P'^^-^y + m]~g{p,m), (16) 



where 

{—i)g{p,m) = (— i) / dsexp[ — trlncosh(eFs) + ij9(eF)~"'^ tanh(eFs)p eF^^a^^s 

Jo 2 2 

p2 _ ^2 _|_ 2 (p2 _ _|_ 4 (^p2 _ ^2 _j_ ^g^3 ^^2 _ ^2 _|_ ^g'j4 

(17) 

Here ^I'^^'^^-'^^'^l is the completely antisymmetrized product (for d = 4, ^[^^'^^•^^'^l = 
To obtain the second line in (IT7|) we have used the identity, |{cr^jy, o"ak} = 
VfiXVyn — Vpti^VyX — Ii^jHuIxIk]- Because of the relation $(x, x')$(x', x) = 1, the phase fac- 
tor $(x, x') is always canceled in charged fermion loops. So the only difference between 
the ordinary Feynman rules and those in an uniform external electromagnetic field is the 
propagator. We need to use the propagator {j) + m)g{p, m) instead of i{f + m)/{p'^ — rn?). 



iiw? — ie)s\ 



Fig. [2] is the lowest order self energy diagram which depends on F^^, 

/d'^k 
^^Tr[(^ + m2)g{k, m2){j) + ^ + mi)~g{j) + k, mi)], (18) 



where Tr means the trace of 7 matrices. We substitute the expansion of g{p, m) in ([Tj 
into (dHD and calculate n(p, F) to 0{F^). As we argued just above eq. ([1]), we consider 
only the terms which depend on F^y because the term which is independent of F^y does not 
contribute to the shift of the pole of the propagator at = /i^. We can calculate n(p, F) 
by the usual techniques for loop calculations [5]. The term linear in F^y vanishes because 
Lorentz invariants F'^ and F^'^p^py are zero. Because the term quadratic in F^y is finite at 
c/ = 4, we will henceforth assume that the regularization has been removed and d = A. To 
calculate the term quadratic in F^j,, we use the identity 

Tr[7^F ■ a-i^F ■ a] = 2^{FY'' - 2V^^^ (19) 

where F ■ a = F^'^^aap and / = —^F^^'^F^y. After some calculation, we obtain the term 
quadratic in Ffj_y 

n(p,F)|^. = + {F'rP,PuB2{p% (20) 

where 

/■I 1 1 

Bi= dx\— {9x{l - x) -2) + —{-mix - ml{l - x) + mim2)], (21) 
Jo A A 

/■I 1 2 

B2= / dx[—{l-x)x{10x{l-x)-l)-—{l-xfx'^{-mlx-ml{l-x)+mim2)],{22) 
Jo ^ 

and A = p'^x{x — l)+mlx + 'ml{l — x). By comparing ([3]) and (l20l) . we obtain the coefficients 
in (ED; a2 = b2 = 0, a, = h = ^|i?i(p' = 0), and 

ci = f^Y^^Mp' = 0). (23) 
We can calculate the x integral in B2{p'^ = 0) and obtain 

mtB2{p^ = 0) = i-^TyMv - i)(z/' + lOy + 1) - Qy{y + 1) Iny] 

(24) 

+ (^3Tj^f^^^ - 1) - (2/ + 1) In?/] + ^^^[3(y2 _ 1) _ (^2 + 4^ + 1) in^], 

where y = m\/m\. In particular, we notice that m2-B2(p^ = 0,?/ = 1) = 1/6 and m\B2{p'^ = 
0, y — )■ 00) ~ — (ln?/)/y^ (we assume mi > m2). So ci in ( 123|1 as a function of y is negative 



m2''B2(p^=0) 
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FIG. 3. mlB2{p^ = 0) in (IMl) as a function oi y = mf/m|. ci = ^^^B2{p^ = 0) can be both 
positive and negative when y varies. 

for large y. We show the graph of m\B2{p^ = 0,y) as a function of y in Fig. [31 and 
we see that Ci can be both positive and negative when y varies. From ([H]) and fllUp . the 
maximum velocities of the scalar particle become greater than unity when ci is negative, so 
the scalar particle in this example can propagate at speed greater than unity in an uniform 
electromagnetic field. 

III. CONCLUSION AND DISCUSSION 

In this paper we studied the maximum velocity of a neutral real scalar field interacting 
with charged fields in an uniform electromagnetic field. We showed that the kinetic terms 
of the scalar fields are modified by the radiative corrections and the maximum velocity of 
the scalar particle becomes greater or less than unity depending on the structure of the 
quantum field theory. From ([8]) and (fTOj) . the maximum velocity become greater (less) than 
unity when ci in ([3]) is negative (positive). The physical interpretation of this effect is that 
the vacuum is considered as the medium filled with imaginary charged particles and the 
external electromagnetic field changes the property of the vacuum. The maximum velocity 
depends on the frame of reference (i.e. the components of F^i,) and causality is not violated 
even if the maximum velocity is greater than unity. From ([HD this change of the kinetic term 
can be considered as the change of the spacetime metric for the scalar field. So the change 
of the kinetic term also modifies probably the lifetime of the particle in an electromagnetic 



We considered the neutral real scalar field in an uniform electromagnetic field to make 



field. 
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the effect of radiative corrections clear. We can generalize the analysis in various ways. For 
example, the maximum velocity of charged fields and of higher spin fields becomes probably 
greater or less than unity depending on the structure of the quantum field theory in an 
uniform electromagnetic field. 

Finally we discuss the possibility of this superluminal propagation in the Standard Model. 
The Yukawa coupling of mesons is greater than one and the perturbative calculation is not 
reliable. So we cannot know the sign of Ci, however, by dimensional analysis we can estimate 
the magnitude of ci; |ci| ~ e^m~^, where rup is the proton mass. From ([8]) the shift of the 
maximum velocity from unity in an uniform electromagnetic field that E and B are parallel 
to each other is \v± — 1\ ^ e^m~'^{Q'^ + E^). As an example, we estimate this value for the 
electric field E„ near a nucleus (e|E„| ^ e^m^): \v_i_ — 1\ ^ ^ 10""^. We can possibly 
observe the change of the maximum velocity in the propagation of light neutral mesons in 
matter. 

A simple case in which the perturbative calculation is reliable is the case of neutrinos . 
The propagator of neutrinos depends on the electromagnetic field through W boson charged 
lepton loops. The mechanism studied in this paper can possibly explain the experiment by 
MINOS collaboration 3] which reported that the speed of neutrinos is greater than unity. 

We have not had an intuitive explanation for the negative sign of ci. For example, we 
can see easily that the sign of Ci for a neutral real scalar field interacting with two charged 
scalar fields which have different masses is always positive. So the difference of masses of 
charged fields is not a sufficient condition. We need extensive studies to know whether there 
is the superluminal propagation in the Standard Model or not. 
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